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The Fermi gas model, while providing a reasonable qualitative description of the continuum nuclear response, 
does not include the effects of dynamical nucleon-nucleon correlations in the initial and final states, that have 
long been recognized to play a critical role in specific kinematical regions. We review a many-body approach in 
which these effects are consistently taken into account and discuss the results of a calculation of the quasielastic 
neutrino-oxygen cross section as an illustrative example. 



1. Introduction 

Within the Fermi gas (FG) model the nucleus is 
seen as a degenerate gas of protons and neutrons 
at density p, binding effects being taken into ac- 
count through the inclusion of an average inter- 
action energy e . 

In spite of its simplicity, the FG model pro- 
vides a remarkably good description of the in- 
clusive electron-nucleus cross section in the re- 
gion of the quasi free peak, corresponding to 
Xb = Q'^ jlmv ~ 1, where xb is the Bjorken scal- 
ing variable, — — v^ , q and v being the mo- 
mentum and energy transfer, respectively, and m 
denotes the nucleon mass. The successful appli- 
cation to the analysis of electron scattering data 
at moderate momentum transfer (q 150 — 200 
MeV) 2 also prompted the extension of the FG 
approach to the case of neutrino-nucleus scatter- 
ing 121 . 

Over the past two decades the limitations of the 
FG description have been exposed by the grow- 
ing availability of new electron scattering data, 
extending into the region of large q and low v, 
corresponding to xb ^- The failure of the FG 
model to explain the measured cross sections in 
this kinematical regime, clearly illustrated in Fig. 
n has been ascribed to the dominance of scatter- 
ing off nucleons belonging to strongly correlated 
pairs, that cannot be accounted for by any mean 
field approach. 

In addition, it has to be pointed out that, as 



the FG model does not include surface and shell 
effects, it cannot be employed to study exclusive 
channels, which are known to be sensitive to the 
details of the nuclear wave function. 
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Figure 1. Cross section versus energy loss for scat- 
tering of 4 GeV electrons at 15 ° from Carbon. The 
dot-dash, dashed and full lines show the quasielas- 
tic, inelastic and total cross sections predicted by the 
PC model with fc^ = 221 MeV and e = 35 MeV. The 
arrow points to the value of v corresponding to quasi- 
free kinematics. The experimental data are taken 
from ref. 2|. 



In this paper we review a theoretical approach, 
based on nonrelativistic nuclear many-body the- 
ory, that incorporates in a consistent fashion 
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dynamical nucleon-nucleon (NN) correlations in 
both the initial and final state. The impulse ap- 
proximation (lA) scheme, which allows one to re- 
late the cross section to the nuclear spectral func- 
tion P{p,E), describing the energy-momentum 
distributions of nucleons in the target nucleus, 
is outlined in Section 2. The main features of 
P{p,E), extracted from both experiments and 
theoretical calculations, are described in Section 
3, while Section 4 is devoted to the discussion of 
final state interactions (FSI). Finally, in Section 
5 we summarize the results and state the conclu- 
sions. 

2. The impulse approximation 

The main assumption underlying I A is that, as 
the space resolution of a probe delivering momen- 
tum q to a nucleus is ~ l/|q|, at large enough |q| 
the target is seen by the probe as a collection 
of individual nucleons. Within this picture, the 
nuclear response measures the probability that, 
after giving one nucleon a momentum q at time 
t = 0, the system be reverted to the ground state 
giving the same particle a momentum — q after 
time t. 

The second assumption involved in lA is that 
(FSI) taking place a,t < t' < t between the hit 
nucleon and the (A-l)-nucleon spectator system 
be negligible. 

In the lA regime the scattering process off a 
nuclear target reduces to the incoherent sum of 
elementary processes involving only one nucleon, 
as schematically illustrated in Fig. |21 



Let us consider the weak charged current reac- 
tion 



iy^ + A^H-+p+{A-l) 



(1) 



According to lA, the corresponding differential 
cross section can be written 



j (fp dE P{p, E) 



da 



N 



dflf^dEf^ 



'.S{E, - E,, - E - Ep>) , (2) 



where daN /dfl^dE^^ is the cross section describ- 
ing the elementary scattering process 



Vy,{k) +n{p) ^ ^ {k')+p{p) 



(3) 



involving a hound nucleon carrying momentum p. 
The spectral function P(p, E), yielding the prob- 
ability of removing a nucleon with momentum p 
from the nuclear ground state leaving the resid- 
ual system with excitation energy E', will be dis- 
cussed in Section 3. 

Up to a kinematical factor, the elementary 
cross section reads 



da 



N 



dQfj^dE^ 



^^L,Ak,k')W>^^{p,q) , 



(4) 



where G = GpcosOc, Gp and 9c being Fermi's 
coupling constant and Cabibbo's angle, respec- 
tively. The lepton tensor L^j, is totally specified 
by kinematical variables, whereas the definition 
of W^'^{p,q), describing the weak interactions of 
a free nucleon, involves the Dirac, Pauli and axial 
form factors Fi, F2 and Fa- 




Figure 2. Pictorial representation of the lA scheme, 
in which the nuclear cross section is replaced by the 
incoherent sum of cross sections describing scattering 
off individual nucleons, the recoiling (A- 1) -particle 
system acting as a spectator. 




Figure 3. Schematic representation of the lA treat- 
ment of the lepton nucleon vertex. The elementary 
scattering cross section is assumed to be the same as 
in free space, binding effects being taken into account 
replacing the energy transfer u with f (see text). 
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Within the lA scheme binding is taken into 
account replacing the physical four-momentum 
transfer q — k^k' = (i^, q) with q= q), where 



M 



A-l 



V\p\ 



(5) 



Ma and Ma~i — Ma — m + E being the tar- 
get mass and the mass of the recoiling nucleus, 
respectively 5 . This essentially amounts to as- 
suming that: i) a fraction (ly — v)/^ of the lepton 
energy loss is spent to put the struck nuclcon on 
the mass shell and ii) the elementary scattering 
process can be described as if it took place in free 
space with energy transfer v < u. 

3. The nuclear spectral function 

Let us consider a system of A nucleons whose 
dynamics is described by the nonrelativistic 
hamiltonian 



i—1 jyi—l 



(6) 



where is the momentum of the i-th nucleon, 
Vij is the NN potential and the ellipsis indicates 
the presence of three-nucleon interactions ^ . 

The spectral function is defined as (see, e.g., 
ref. El) 

Pip.E) = E|(*n^''Vpl*^)f 

n 

xSiE + Eo-En), (7) 

where l^*^) is the nuclear ground state, satisfy- 
ing the Scrodinger equation HaI'^q) — Eo\'^^), 
while and En denote the n-th eigenstate 

of the (A-l)-nucleon system and the associated 
energy eigenvalue, respectively. 

Due to the complex spectrum of (A-l)-nucleon 
states, in general the calculation of P(p, E) 
within nuclear many-body theory involves pro- 
hibitive difficulties. It has been carried out only 
for few- nucleon systems, having A= 3 |7I8| and 
4 [3, and uniform symmetric nuclear matter, i.e. 
in the hmit A ^ oo with Z=A/2 pTHj . 

'^The inclusion of a three-nucleon potential in nuclear 
many-body theory is needed to reproduce the binding en- 
ergy of the three-nucleon system and the empirical equi- 
librium density of nuclear matter. 



Spectral functions of nuclei with A > 4 have 
been constructed using the Local Density Ap- 
proximation (LDA) JT]. Within this approach 
P(p, E) is split into two parts, corresponding to 
low momentum nucleons, sitting in single parti- 
cle shell model states, and high-momentum nu- 
cleons, belonging to strongly correlated pairs, re- 
spectively. 

The spectral function of the shell model states 
has been extensively measured by single nu- 
cleon knock-out experiments (see, e.g., ref. jOj), 
whereas the correlation part is obtained from the 
results of theoretical calculations through 

P^orrip, E)= I d\ p{v)PZ^!.{p, E; p{r)) (8) 

where p{r) is the nuclear density distribution and 
^ra*/(Pi -E-; p) is the correlation part of the spec- 
tral function of infinite nuclear matter at uniform 
density p. 
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Figure 4. Three-dimensional plot (left panel) and 
scatter plot (right panel) of the oxygen spectral func- 
tion obtained using the LDA approximation described 
in the text. 



The LDA spectral function of ^^O obtained 
combining the nuclear matter results of ref. [HI 



4 



and the Saclay (e, e'p) data TT is shown in Fig. 
^ The sheU model states account for 80 % of 
the strength, whereas the remaining 20 % is 
located at high momentum (|p| 3> kp) and large 
removal energy {E ^ Ep, where Ep denotes the 
Fermi energy) . It has to be emphasized that large 
E and large p are strongly correlated. For exam- 
ple, ^ 50 % of the strength at |p| = 320 MeV is 
located at > 80 MeV. 

FigiElshows the momentum distribution of nu- 
cleons in the ground state of ^^O, obtained from 
the LDA spectral function through 



n(p)= J dE P{p,E) 



(9) 



Comparison with the results of a Monte Carlo 
calculation carried out using a highly realistic 
many-body wave function 14 suggests that LDA 
results are quite reasonable. For reference, the 
Fermi gas momemtum distribution corresponding 
to kp = 220 MeV is also shown by the dashed line. 



-LDA 

<J Monte Carlo calculation 

(S.C. Pieper) 
-FG ( kp = 221 MeV ) 




Figure 5. Momentum distribution of nucleons in the 
oxygen ground state. Solid line; LDA approximation. 
Dashed line: Fermi gas model with kp — 220 MeV. 
Diamonds: Monte Carlo calculation carried out by 
S.C. Pieper using the wave function of ref. |14| . 



pi~ + p + X is displayed in Fig. for the case of 
neutrino energy Ei^ = 2 GeV and muon scatter- 
ing angle 30°. Comparison between the solid line, 
obtained using the LDA spectral function, and 
the dashed line, corresponding to the FG model, 
shows that the inclusion of dynamical NN correla- 
tions produces a quenching of the peak associated 
with the appearance of sizeable tails, particularly 
at large lepton energy loss. 
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Figure 6. Cross section of the process +^'^ O 

+ p + V. Solid line: lA calculation carried out 
using the LDA spectral function discussed in the text. 
Dashed line: Fermi gas model with kp = 220 MeV 
and 1=25 MeV. 



The presence of high-momentum and high- 
energy components in the nuclear wave function 
is also responsible for the shift in the position of 
the peaks. The average removal energy 



{E) = J d^p dE EP{p, E) 



(10) 



calculated using the LDA spectral function of 
-^^O, {E) = 42 MeV, is in fact much larger than 
the binding energy Z — 2b MeV of the FG model. 



4. Final state interactions 



The lepton energy loss dependence of the cal- 



culated cross section of the process v^^ 



O 



The occurrence of FSI leads to the breakdown 
of the lA picture of the scattering process, thus 
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making the reconstruction of its kinematics much 
more difficuh. 

A number of theoretical studies of FSI in 
(e, e'p) have been carried out within the Dis- 
torded Wave Impulse Approximation jl5j . in 
which the spectral function is modified to take 
into account the distortion of the struck nucleon 
wave function produced by its interactions with 
the spectators, described by a complex optical po- 
tential. DWIA has been succesfuUy employed to 
analize two-body breakup processes, in which the 
spectator system is left in a bound state, for pro- 
ton energies IS 100 MeV. 

Realistic many-body calculations of nuclear 
structure provide strong a priori indication that 
the PWIA picture is not likely to be applicabile 
even at much lager energies. 
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Figure 7. The function {Npi{h,z)), defined by 
Eg. dllll . plotted at 6 = 0.25 fm for a 1 GeV nucleon 
knocked out from oxygen. 



The probability that a nucleon hit at position 
r in the target center of mass frame and left with 
momentum p' in the direction of the z-axis rescat- 
ter against the spectator particles is illustrated in 
Fig. 13 showing the z-dependence of the quantity 

(A^p'(b,^)> dR |vI/o(i?)P^<5(r-r,) 



ctnn(p') 



case of a 1 GeV nucleon knocked out from oxygen. 
In Eq.(|ni p(j) and *o(^), with i? = (ri, . . . , r^), 
denote the target density and ground state wave 
function, respectively, whereas (TjvAf(p') is the NN 
scattering cross section at incident momentum p'. 
The 3A-dimensional integral has been carried out 
using Monte Carlo configurations sampled from 
the probability distribution associated with the 
^^O wave function of ref. |14| . 
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Figure 8. Schematic representation of rescattering in 
nucleon knock-out processes. 



The quantity (A'p' (b, z)) is the average number 
of spectators found within the cylindrical volume 
described by a circle of area crA'Ar(p') initially cen- 
tered at r and moving along the direction of the 
momentum p' (see Fig. Its value provides 

an estimate of the number of rescatterings un- 
dergone by the knocked out particle. The results 
shown in Fig. [7| suggest that in the case of scat- 
tering processes involving nucleons in the z < 
emisphere, the struck particle is likely to interact 
with at least one of the spectators. 

A theoretical treatment of the corrections to 
the lA cross section due to FSI effects was de- 
veloped in ref. |lfi| and extensively applied to 
the analysis of inclusive electron-nucleus scatter- 
ing data |9I11I16| . The main approximations in- 
volved in the approach of ref. ^S] , expected to be 
reasonable at large p', are that i) the struck parti- 
cle moves along a straight trajectory with a con- 
stant velocity v {eikonal approximation) and ii) 
the spectator particles are seen by the fast struck 
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nucleon as a collection of fixed scattering centers 
(frozen approximation) . 

Under the above assumptions the ground state 
averaged propagator of the struck nucleon can be 
written in the factorized form 



0.5 



C/p.(t) = C/°,(t)[/™(i), 



(12) 



where Up, is the free space propagator, while FSI 
effects are described by the quantity 



i=l 



) -(13) 



In Ea. fT^ Tij = ri — Tj, Wp'(|r|) is the coordinate- 
space NN scattering t-matrix, parametrized in 
terms of total cross section, slope and real to 
imaginary part ratio, and (...) denotes the ex- 
pectation value in the target ground state. 

Note that U^,^^ {t) is simply related to the nu- 
clear transparency Tp/, measured in coincidence 
(e, e'p) experiments QT], through 



Tp,= hm \Up\t)\^ 



(14) 



In presence of FSI the inclusive cross section 
can still be expressed in terms of the lA result of 
Eq.© through [El 



da 



- / dE' 



^ dn^dE'^ 



- U,{E^-E'^) , (15) 



the folding function fp'ito) being defined as 



dt 
2^ 



1/2 



FSI 



it)-T;i 



1/2 



(16) 



The above equations clearly show that the 
strength of FSI is measured by both Tp, and 
the width of the folding function. In absence of 
FSI Up,^^{t) = 1, implying in turn Tp, = 1 and 
/p'(w) -^S{uj). 

1 /2 

Fig. El shows the quantity /p'(w) — S{Lo)Tp', , 
evaluated for a 1 GeV nucleon knocked out from 
oxygen using a parametrization of the measured 
NN scatering amplitude ^Sl and the wave func- 
tion of ref. '14'. 

The effect of FSI on the calculated cross section 
of the process ly^ O ^ fi^ +p + X a.t neutrino 




-0.4 -0.2 0.0 0.2 
cj [GeV] 



0.4 



Figure 9. Finite width contribution to the folding 
function of Eg. (1161 . evaluated for a 1 GeV nucleon 
knocked out from oxygen using a parametrization of 
the measured NN scattering amplitude |18| and the 
wave function of ref. |14|. 



energy E^ — 2 GeV and muon scattering angle 
30° is illustrated in Fig. [TOI It clearly appears 
that, in spite of the fact that the finite width 
component contributes only ^ 15 % of the fold- 
ing function normalization, the inclusion of FSI 
produces a sizable redistribution of the strength, 
resulting in a quenching of the peak and an en- 
hancement of the tails. 

5. Summary and conclusions 

Nuclear many-body theory provides a consis- 
tent framework allowing for a description of the 
continuum nuclear response that takes into ac- 
count dynamical effects not included in the FG 
model. 

The approach discussed in this paper, that 
has been extensively applied to the analysis of 
electron-nucleus scattering data, is based on the 
lA picture, in which the lepton probe scatters 
off individual nucleons distributed in momentum 
and energy according to the spectral function 
P{p,E). 

While this picture proved to work exceedingly 
well for scattering of few GeV electrons, its appli- 
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Figure 10. Cross section of the process 
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/x + p + X. Solid line: Full calculation, including FSI. 
Dashed hne: lA calculation. Dot-dash line: Fermi 
gas model with A;f = 220 MeV and e = 25 MeV. 



cability in the region of lower energies relevant to 
neutrino experiments needs to be quantitatively 
investigated. For example, the effect of Pauli 
blocking of the struck nucleon is likely to play 
a non negligible role at Ey < 700 MeV. 

Comparison between the results of the FG 
model and those obtained using a realistic spec- 
tral function shows that the inclusion of dynami- 
cal correlation effects, leading to the appearance 
of high momentum and high removal energy com- 
ponents, sizably affects the response. 

Corrections to the lA cross sections due to the 
occurrence of FSI, which have long been recog- 
nized to play a major role in redistributing the 
inclusive strength, can also be taken into account 
in a consistent fashion within nuclear many-body 
theory. 

As a final remark, it s worthwhile mentioning 
that the proposed many-body approach can be 
readily generalized to include the contribution of 
inelastic channels in a fully consistent fashion. 
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